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Hydrodynamic expansion of the hot fireball created in relativistic Au-Au collisions at y/s = 
200GeV in 3 + f-dimensions is studied. We obtain a simultaneous, satisfactory description of the 
transverse momentum spectra, elliptic flow and pion correlation radii for different collision cen- 
tralities and different rapidities. Early initial time of the evolution is required to reproduce the 
interferometry data, which provides a strong indication of the early onset of collectivity. We can 
also constraint the shape of the initial energy density in the beam direction, with a relatively high 
initial energy density at the center of the fireball. 
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I. INTRODUCTION 

A multitude of experimental data from the Relativis- 
tic Heavy Ion Collider (RHIC) indicate that dense, col- 
lectively expanding matter is created in ultrarelativis- 
tic nuclear collisions [l|, HI, HI, 0] ■ Ratios of multiplicities 
of different particles produced in central collisions can 
be described assuming chemical equilibration of parti- 
cle abundances [1, HI, [3, HI, HI, Hi] • Transverse momen- 
tum spectra of particles produced at central rapidities 
are thermal up to transverse momenta of about 2GeV/c. 
Particle spectra result from a collective, transverse ex- 
pansion of the matter coupled with subsequent thermal 
emission jTl], EH- Possible rescattering and resonance 
decays have been modelled and the conclusion, that at 
some stage of the expansion a dense locall y eq uilibrated 
fireball is formed, remains unchanged fl3l, Il4|. Another 
measured quantity directly resulting from the collective 
expansion is the elliptic flow. For non-zero impact pa- 
rameters the fireball is azimuthally asymmetric in the 
transverse plane, and its expansion imprints the momen- 
tum distribution of final hadrons with measurable az- 
imuthal asymmetry [ill, [l(| [l?], HH ■ The existence of the 
dense matter is demonstrated in yet another way by the 
observation of the attenuation of the production of high 
energy hadrons in nuclear collisions. This effect is due to 
the energy loss of energetic partons while traversing the 
dense fireball O M, El, HI HI Hi • 

Relativistic hydrodynamics is very well suited for the 
description of the collective phase of the fireball expan- 
sion Q HE Hi, Hi, H3, Hi, Hi, Hi, HH, Hi, Hi, Hi, HI Hi, 

H3- Assuming local thermal equilibration, perfect fluid 
hydrodynamics can be used. Starting from an initial en- 
ergy density profile, the fluid expands and cools down. 
In the process, gradients of the pressure cause the accel- 
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eration of the fluid elements and collective flow velocity 
is formed. In the longitudinal (beam) direction Bjorken 
flow with velocity v z = z/t is usually assumed in the ini- 
tial conditions. On the other hand, the appearance of a 
substantial transverse flow can be considered as a robust 
signature of the formation of strongly interacting mat- 
ter in the overlap region of heavy-ion collisions. Most 
of the hydrodynamic calculations modeling nuclear col- 
lisions at RHIC energies assume boost-invariance [38| in 
the beam direction. Such approaches are effectively 2+1- 
dimensional (2 + ID) and are restricted to central rapidi- 
ties. Existing experimental data outside of the central 
rapidity region on particle multiplicity and elliptic flow 
show that at RHIC energies the Bjorken boost-invariance 
is not realized. Calculations exist for the general 3 + ID 
geometry of the collision 0, H3, Hi, M, H|- They 
show that relativistic hydrodynamics can be applied for 
a broad range of rapidities in central and semiperipheral 
collisions. These studies can describe transverse momen- 
tum spectra and the elliptic flow of produced particles. 
On the other hand, Hanbury Brown-Twiss (HBT) corre- 
lations between identical particles cannot be accounted 
for [27], Hi, Hi, HH • In the present paper we investigate 
3 + ID hydrodynamic expansion of the fireball and show 
that a simultaneous and satisfactory description of the 
particle spectra and elliptic flow (for a broad range of 
rapidities) as well as of the HBT radii for central rapidi- 
ties can be achieved. The key ingredients of the model 
leading to this success are the use of a realistic equation 
of state without a first order phase transition and a rel- 
atively early start up time for the collective expansion. 
This hard equation of state and the small initial time in- 
dicate that the initial state is a highly compressed matter 
with energy density of up to 100GcV/fm 3 in central col- 
lisions. In this work we use perfect fluid hydrodynamics. 
Shear viscosity or hadronic dissipative effects are known 
to modify final observables [H HI Hi, Hi, S, H ES E3 , 
especially the elliptic flow. The influence of viscosity ef- 
fects on particle spectra or HBT radii is more difficult 
to be explicitly demonstrated, since such effects can be 
compensated by a change in the unknown initial time or 
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energy density profile. 



II. HYDRODYNAMIC EQUATIONS AND 
INITIAL CONDITIONS 



In a perfect fluid each element is locally in thermal 
equilibrium. At each point the fluid is characterized by 
its four velocity it' 1 , the energy density e, and the pressure 
p. The energy momentum tensor is 



T" v = (e + p)u^u v 



9^ P 



Hydrodynamic equations 



= 



(2.1) 



(2.2) 



in the full 3+ ID geometry represent 4 independent equa- 
tions, and together with the equation of state allow to 
calculate the evolution of the densities and velocities of 
the fluid starting from some initial conditions. We use 
a realistic equation of state interpolating between lat- 
tice data at high temperature (above the critical tem- 
perature T c = 170MeV) and an equation of state of a 
nonintcracting gas of massive hadrons at lower temper- 
atures 48]. This equation of state presents only a very 
moderate softening around the critical point. The use 
of this realistic equation of state is the key to the suc- 
cess of 2 + ID hydrodynamic description of RHIC data 
on transverse momentum spectra, elliptic flow and HBT 
radii [2(| [4f| . For the modelling of the expansion of the 
fireball created in ultrarelativistic collisions it is useful to 
define the proper time and space-time rapidity variables 



t + z 
t- z 



(2.3) 



with z the beam axis coordinate. The four velocity is pa- 
rameterized using the two components of the transverse 
velocity u x and u y and the longitudinal fluid rapidity Y 



= (^fcoshY,u x ,Uy,'ysmhY) 



(2.4) 



where 7 = « / 1 + u 2 x + . Hydrodynamic equations re- 
late four unknown functions : the velocity fields Y , u x , 
and u y and either the energy or the pressure. In practice 
the numerical solution is more stable if instead of the en- 
ergy density (pressure) the logarithm of the temperature 
is used 



T = \og{T/T L ) 



(2.5) 



with Tjr a constant temperature. The velocities and T 
are function of r, 77, x, y and the hydrodynamic equations 



can be written in the following form 



DT 



7 ( sinh(r - V )dr + cosh(y V \ ) Y 



uv s A . sinh(Y-77) 
cosh(r — rj)o T H a v 7 



+ X U X + dyUy 

Du x = -(1 
—v-xl ( cosh(F — rj)d T 



U 2 x )d X T - U X Uy3yT 

sinh(y — if) 



Du,. 



-(1 + u 2 y)d y T - u x u y d x T 



u y j cosh(F — T])d T + 



D(7sinhy) = — coshF 



sinh(Y — rj) 



0, 



sinh(y — rj)d T 
sinh(y — 



cosh(F — 77) 



+ (u\ + u 2 y )(cosh(Y - T])d T 

— 7sinhF (u x d x + u y d y ) T , 
where c s is the sound velocity and 



T 



(2.6) 



D 



u^d^ = u x d x + Uydy 



+ 7(cosh(F — rf)d T + 



sinh(y — 77) 



d v ). (2.7) 



The first of Eqs (|2.6[) is the entropy conservation equation 
d„{u»s) = 0. 

The differential equations (|2.6[) are solved as an evolu- 
tion in proper time starting from some initial conditions 
at r = to. At the initial time there is no transverse flow 
(u x = and u v — 0), the initial longitudinal rapidity 
follows the Bjorken scaling flow Y(tq, 77, x, y) — 77. Early 
initial time of the hydrodynamic evolution tq = 0.25fm/c 
implies a high energy density in the initial state. The sys- 
tem evolves for a longer time, which leads to a stronger 
transverse as well as longitudinal flow. Experimental ob- 
servation of a strong rapidity dependence of the elliptic 
flow and of the particle densities [5l[ indicates that 
the Bjorken scaling scenario is not realized at RHIC en- 
ergies. It means that a Bjorken scaling plateau in the 
initial energy density distribution in space-time rapidity 
cannot extend over a large interval. In the transverse 
plane the energy density is assumed to be proportional 
to a combination of Glauber Model densities of wounded 
nucleons and binary collisions. The initial energy density 
distribution at impact parameter b is 

e(ro) - kf(ri-r) sh )[{N A (x,y)+N B (x,y))(l-a) 

+ 2aN bin (x,y)] . (2.8) 

Na and Nb are the densities of wounded nucleons from 
the right and left moving nuclei respectively, N^n is the 
density of binary collisions 

N A {x,y) = T(x-b/2,y)(l-exp(-aT(x + b/2,y)/A)) 
N B (x,y) - T(x + b/2,y)(l-exp(-<jT(x-b/2,y)/A)) 
N Hn {x,y) = aT(x-b/2,y)T(x + b/2,y) (2.9) 



3 



and 



T ( x , V) = J dzp(x, y, z) (2.10) 



is the thickness function calculated from the Woods- 
Saxon density of colliding nuclei 



Pa 



1 + cxp {[[^/x 2 + y 2 + z 2 - R A )/a 



(2.11) 

For Au nuclei (A = 197) we take p = 0.1 7fm 6 , R A = 
6.38fm and a — 0.535fm; the inelastic cross section is 
a = 42mb. The density of wounded nucleons N A + N B 
is used to calculate the total number of participants at 
each impact parameter, these numbers are used to fix the 
impact parameters corresponding to centrality bins used 
in the analysis of the experimental data. The profile in 
the longitudinal direction is 



f(v) = ex P 



2 tf 



\m -vo. 



(2.12) 



with a plateau of width 2rjo = 2.0 units in space-time 
rapidity, and Gaussian tails with half width a v = 1.3. 
At each point in the transverse plane the distribution 
in space-time rapidity is shifted by the center of mass 
rapidity of the local fluid [27| 



Vsh = ^ log 



N a + Nb + vn(N a -N b ) 
N A + N B -v N {N A -N B ) 



(2.13) 



where vn is the velocity of the projectile in the center of 
mass frame. The coefficient of k in Eq. (|2.8[) is taken so 
that the energy density at the center of the fireball at zero 
impact parameter is 107GeV/fm 3 at To = 0.25fm/c. It 
corresponds to a temperature of 510MeV, well above the 
critical temperature. The initial distributions for other 
impact parameters are obtained from geometrical scaling 
(|2.8p only, with a contribution of binary collisions a — 
0.145 [52j. This provides a satisfactory description of 
charged particle multiplicities for centralitics — 40%. 



III. EVOLUTION OF THE HOT MATTER 



The numerical solution of Eqs. (|2.6[) is obtained as an 
evolution in proper time from initial densities (|2.8p . The 
total entropy (s is the entropy density) 



S = J 7Cosh(y — r/)s dxdxdrj 



(3.1) 



is conserved to the accuracy of less than 0.5%. At the 
very beginning of the evolution a very rapid longitudinal 
expansion occurs (Fig. [1]) . The matter at the edges of the 
plateau of the energy density distribution is subject to 
longitudinal acceleration and eventually the distribution 
becomes approximately a Gaussian, that grows wider in 




-6 -4 -2 



FIG. 1: (Color online) Energy density as function of space- 
time rapidity for different proper times r = 0.25, 2, 4, 6fm/c 
(dashed-doted, solid, dotted and dashed lines respectively). 
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FIG. 2: (Color online) Deviation of the longitudinal fluid 
rapidity from the Bjorken flow Y(r, rj, x, y) — r\ for different 
proper times r = 2, 4, 6fm/c ( solid, dotted and dashed lines 
respectively) . 



time. The longitudinal acceleration is known to depend 
on the equation of state [53| and on possible viscosity 
effects [54]. For the perfect fluid and a hard equation 
of state the longitudinal expansion and acceleration is 
significant [53]. 

For comparison we calculate a hydrodynamic evolu- 




5.0 10.0 

t [fm/c] 

FIG. 3: (Color online) Temperature at the center of the fire- 
ball as function of the proper time from the 3 + ID (dashed- 
dotted line) and from the 2 + ID (solid line) evolutions 
(b = 2.1fm). 
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FIG. 4: (Color online) Freeze-out hypersurface Tf = 150MeV, 
for the impact parameter b — 2.1fm in the plane (t — x), y = 0, 
for rj = (dashed line) and r\ = 3 (dotted line). The solid line 
represents the freeze-out hypersurface for the 2+1D evolution. 
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time rapidity in the initial conditions ()2 . 13[) . In Figs [4] 
and [5] is shown a cut (r — x) through the freeze-out hy- 
persurface at the freeze-out temperature Tf = 150MeV. 
At the central space-time rapidity 77 = the freeze-out 
hypersurfaces in the 3 + ID and 2 + ID calculations are 
very similar. For central collisions the dense system ex- 
ists for lOfm/c. This short life-time of the fireball results 
in values of extracted HBT radii compatible with the ex- 
periment. For large space-time rapidities the transverse 
size and the life-time of the fireball is smaller. The asym- 
metric shape of the freeze-out hypersurface for 77 7^ is 
the result of the tilt in the initial flow of the matter (Eq. 

Most general freeze-out hypersurfaces realized in the 
hydrodynamic expansion can be parameterize using 3 an- 
gles 



ths = = 

Vhs = 

xhs = 
Vhs = 



d(9, £, 4>) sin £ sin 9 + tq 



A 



• cos 9 



d(9, £, 4>) cos £ sin 9 cos 4> 
d(9, £, <fi) cos £ sin 6 sin t/> ; 
< 9 < 7T 
< £ < tt/2 
< <p < 2tt , 



(3.2) 



FIG. 5: (Color online) Same as Fig. g]but for b = 7.1fm. 



tion of the system assuming a 2 + ID boost invariant 
expansion, with the energy density profile in the trans- 
verse plane given by equation (|2.8[) . but without 77 de- 
pendence. The temperature at the center of the fireball 
(T = 485MeV for 6 = 0) that reproduces the observed 
spectra is slightly lower than in the 3 + ID case. This 
is the effect of the additional cooling in the non-boost 
invariant geometry due to the longitudinal acceleration. 
In 3 + ID the longitudinal fluid rapidity is larger than 
in the Bjorken scaling solution (Fig. [2]). At the center 
of the fireball the temperature drops down following the 
Bjorken formula T oc r~ c = (Fig. [3]) up to r = 2 — 3fm/c. 
Later cooling from the transverse expansion and in the 
case of 3 + ID additional longitudinal colling set in. As 
a result the life-time of the 2 + ID and 3 + ID systems 
is very similar, in spite of the fact that in the later case 
the initial energy density is a factor 1.25 higher. 

The hydrodynamic evolution is followed until freeze- 
out, that is assumed to happen at fixed temperature, 
with particles emitted from the freeze-out hypersurface 
without further rescattering. In the following, we present 
results for two different freeze-out temperatures Tf = 150 
and 165MeV. The freeze-out hypersurface is a three- 
dimensional surface in T,T],x,y coordinates. Its shape 
can be deformed due to a strong collective flow and is de- 
formed in the 77 — a; direction due to the shift in the space- 



A is a constant length. 

Following the Cooper- Frye prescription [55| . particle 
spectra are given by 



d 3 N r 

B^j = J dS^/(PM«") • (3-3) 



dSfj, = e ^vapdex" d^x" dfjtXp d9d(,d(j) is the integration el- 
ement on the freeze-out hypersurface and / is the equi- 
librium Bose or Fermi momentum distribution. The four 
momentum of the emitted particle is 



= {m±y,p±cos<pp,p±sm(pp,m±y) , (3.4) 



and 



p^u^ = mj_7cosh(y — y) — p±(u x cos<p p + u y sin<^ p ) . 

(3.5) 
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FIG. 6: (Color online) Pseudorapidity distribution of charged 
particles for centrality classes — 6%, 6 — 15%, 15 — 25%, 
25 - 35%, 35 - 45% and 45 - 50% calculated for the freeze- 
out temperatures Tf = 165 and 150MeV (solid and dashed 
lines respectively) compared to PHOBOS Collab. data (dots) 
[57j . The squares rep resent the BRAHMS Collab. data for 
centrality - 5% [nf. 
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After the hydrodynamic evolution, the 3 dimensional hy- 
persurface parameterized by the variables 9,£,4> is ex- 
ported to the statistical emission and resonance decay 
code THERMINATOR [56]. The density (O (with 
d3HD and (jniSD) is implemented in the code. THERMI- 
NATOR generates events in two steps. First 380 different 
kind of particles and resonances emitted from the hyper- 
surface are generated according to the density Q3.3p . then 
resonances are allowed to decay. 



IV. PARTICLE SPECTRA, FLOW, 
CORRELATION RADII 



Charged particle 
VPS = 



distributions in pseudorapidity 
have been measured for different 



FIG. 7: (Color online) Transverse momentum spectra of tt + 
^ for different rapidity windows ranging from —0.1 < y < 
to 3.4 < y < 3.66 for centrality — 5% (results for different 
rapidity bins are scaled down by powers of 1/10). The dots 
represent the data of the BRAHMS Collab. 

)sinC) 

centralities. The 3 + ID hydrodynamic model can re- 
produce the data for centralities — 40% (Fig [5]). We 
HS show results for two freeze-out temperatures Tf = 165 
and 150MeV. The first one is close to the estimate of 
the chemical freeze-out temperature in Au-Au collisions 
[E S 111- When decreasing the freeze-out temperature 
particle multiplicity goes down, but the effect is small. 
It gives confidence to our model, that assumes a chemi- 
( 3 -£ally equilibrated fluid down to Tf = 150MeV. For lower 
freeze-out temperatures the difference between a chem- 
ically equilibrated and a partially equilibrated fluid be- 
comes significant [2?} ■ The centrality dependence that 
we predict comes solely from the geometrical scaling of 
the fireball density according to Eq. (|2.8[) . other param- 
eters (in particular the freeze-out temperature) remain 
unchanged. On general grounds, one expects the hydro- 
dynamic model to break down for very peripheral col- 
lisions. The interaction region in peripheral collisions 
is not dense enough to equilibrate completely. Experi- 
mental data on the centrality dependence of strangeness 
production and of particle spectra suggest that at impact 
parameters b > 9fm [59, 60, 6l[ less than 70% of the inter- 
action region can be treated as a thermally equilibrated 
fireball. 

For central collisions c = — 5% we calculate the trans- 
verse momentum spectra of pions and kaons at differ- 
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FIG. 8: (Color online) Transverse momentum spectra of K + 
for different rapidity windows ranging from —0.1 < y < 
to 3.2 < y < 3.4 for centrality — 5% (results for different 
rapidity bins are scaled down by powers of 1/10). The dots 
represent the data of the BRAHMS Collab. [|l]|. 
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FIG. 9: (Color online) Transverse momentum spectra of ix + 
for different centrality classes — 5%, 5 — 10%, 10 — 15%, 
15-20%, 20-30%, 30-40% and 40-50% (results for different 
centralities are scaled down by powers of 1/10) calculated for 
the freeze-out temperatures Tf = 165 and 150MeV (solid and 
dashed lines respectively) compared to PHENIX Collab. data 
(dots) [scl- 



ent rapidities. In the whole rapidity range where iden- 
tified particle spectra are available [5l|, we find an ex- 
cellent agreement between the results of the hydrody- 
namic evolution coupled with statistical emission and the 
BRAHMS Collab. data (Figs. Eland©. The best agree- 
ment is achieved for a freeze-out temperature of 150MeV. 
This very good agreement between the data and the hy- 
drodynamic model indicates that the matter created in 
central collisions behaves as a thermally equilibrated (in 
the rapidity range —3.5 < y < 3.5) although not boost- 
invariant fireball. 

A different way of testing thcrmalization in Au- Au col- 
lisions is to compare predictions and experimental data 
for transverse momentum spectra at different centrali- 
ties. In Fig. [9] are shown ir + spectra for p± up to 
3GeV/c and centralities in the range — 50%. We find 
that hydrodynamic calculations with Tf — 150MeV are 
in very good agreement with the experiment for all cen- 
tralities and transverse momenta p± < 2GeV/c. Pions 
with higher transverse momenta originate mostly from 
hard processes and cannot be described as emitted ther- 
mally from a collectively expanding fluid. The agreement 
between the experimental spectra and the results of the 
calculation for K + (Fig. ITU]) is limited to centralities 
- 30%. For centrality bins 30 - 40% and 40 - 50% 
the calculation overpredicts the kaon multiplicity, even 
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FIG. 10: (Color online) Same as Fig |9] but for K+ . 
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FIG. 11: (Color online) Pseudorapidity dependence of the 
elliptic flow coefficient for charged particles for centralities 
15 — 25% for freeze-out temperatures Tf = 150 and 165MeV 
(dashed and solid lines respectively), data for the PHOBOS 
Collab. are denoted by dots [Hoi ]. 
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FIG. 12: (Color online) Transverse momentum dependence of 
the elliptic flow coefficient for protons (dotted line and dots) 
and for 7r + and K + (solid line and triangles). Calculations 
are performed for Tf = 150MeV, data are from the PHENIX 
Collab. HI. 
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FIG. 13: (Color online) HBT radii for Au-Au collisions at 
centrality - 5%. 3 + ID calculations with Tf = 165MeV 
(solid line) and Tf = 150MeV (dashed line), 2+ ID calculation 
with T f = 150MeV (dotted line) and STAR Collab. data [HI 
(squares) are shown. 



though the slope of the spectra is similar in the model 
and in the data. The reduced strangeness production in 
peripheral collisions can be an effect of energy and mo- 
mentum conservation [63| . canonical suppression ImL o r 
reduced size of the thermally equilibrated fireball |59l.l6l|. 
In Figs. 151 and [POl are also shown the results of a 2 + ID 
hydrodynamic calculation for two centralities — 5% and 
20 — 30% (dotted lines, indistinguishable from the 3 + ID 
results). The resulting spectra for central rapidities are 
very similar to the ones from the 3 + ID calculations, 
but are obtained after the expansion of the matter with 
smaller initial energy density. 

The elliptic flow represents a very sensitive probe of 
the collective behavior of the dense matter [18j. The 
azimuthal asymmetry with respect to the reaction plane 
is described by the elliptic flow coefficient v-i 



dN 



dN 



p±dp±d(f) p p±dp_ 



■(l + 2u 2 (p x )cos(2<^) + ...) • (4.1) 



The elliptic flow coefficient for charged particles has 
been measured for a broad range of pseudorapidities [50| , 
showing a strong pseudorapidity dependence. There is no 



sign of a Bjorken plateau for central rapidities. To repro- 
duce the shape of the vi pseudorapidity dependence, ini- 
tial conditions in energy density with a relatively narrow 
plateau in space-time rapidity must be chosen (Eq. I2.12p 
[l4l . [33l ]. Such initial conditions, combined with a hard 
equation of state and an early initial time of the evolution 
result in a complete disappearance of the Bjorken plateau 
in the final hadron distributions. It must be stressed how- 
ever, that it does not mean that for non-central rapidities 
the evolution is not described by hydrodynamics and sta- 
tistical emission. The model model works very well and 
describes the observed spectra for —3.5 < y < 3.5 (Figs. 
[7] and [8]). The longitudinal expansion and the smaller 
size of the system at non-zero space-time rapidities re- 
duce the final elliptic flow. The elliptic flow as function 
of p± for identified particles is shown in Fig. 1121 Hydro- 
dynamic calculations describe the elliptic flow for mesons 
with p± < 1.5GeV/c, but to reproduce the saturation of 
t>2 for large p± dissipative or viscosity effects must be in- 
voked. The elliptic flow for baryons is overpredicted for 
the freeze-out temperatures chosen. 

Pairs of identical particles emitted from the thermal 
source can be used to extract the size of the fireball from 
the interferometry measurement [6?], [H, [69|, [7(|. The 
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statistical emission code THERMINATOR [II] provides 
the space-time points of particle creation. For a set of 
generated events the correlation function is constructed 
from same-event and mixed-event pairs (7l| . Such a gen- 
eral procedure allows for an easy implementation of ex- 
perimental cuts, final interaction or Coulomb corrections 
|72[ [73| and can be also applied to non-identical particle 
correlations [74| . The extracted multidimensional corre- 
lation functions are parameterized by the Bertsch-Pratt 
formula [zllzil 

C(k±,q out ,q side ,qi ong ) = 1 + A exp(- R^ ide q^ ide 

The fit parameters R (HBT radii) are extracted for fixed 
bins of pair total transverse momenta k±. The experi- 
mental radii R out , R S i de and Ri ong are well described by 
the hydrodynamic calculations with Tf = 150MeV (Fig. 
IT3]) . We notice that the 3 + ID and 2 + ID calculations 
that both describe the transverse momentum spectra, 
give also very similar HBT radii. Similarity in the spectra 
means that the transverse collective flow is the same; to- 
gether with the similarity in the freeze-out hypersurfaces 
(Fig. U]) it explains why the HBT radii from the 3 + ID 
and 2 + ID evolutions come out so close. The discrep- 
ancies between the calculations and the data are smaller 
than 10%. The ratio Rout/Rside from the model comes 
out close to the data as well. This remarkable property of 
modern hydrodynamic calculations, which solves the so 
called RHIC HBT puzzle, has been noticed in Ref. (2H ] 
based on 2 + ID simulations. The present study is the 
first 3 + ID calculation using the same equation of state 
and an early initial time. 

V. CONCLUSIONS 

We present an extensive study of the 3 + ID hydrody- 
namic model of the evolution of the fireball. Compared 
to other similar calculations we use a different equation of 
state [48j and an early initial time to start up the expan- 
sion. As a result we find a very good agreement between 
the calculated transverse momentum spectra for pions 
and kaons for different centralities and a broad range of 
rapidities. The fireball formed in Au-Au collisions at 
y 7 ^ = 200GeV is thermalized for centralities — 40%. 
The result is remarkable, since we use an energy density 
profile for different impact parameters predicted from the 



Glauber Model geometrical scaling, with an overall nor- 
malization fixed for most central collisions. It is not sur- 
prisin g th at for peripheral collisions this scaling breaks 
down [60|| , and we overpredict the size of the thermal- 
ized fireball. The thermal description of the production 
of strange particles is subject to even more restrictions. 
As a result the spectra of kaons are well reproduced for 
centralities — 30%. For more peripheral collisions the 
number of observed kaons is smaller than predicted in the 
model, signaling an incomplete chemical equilibration of 
the interaction region. 

The caculation reproduces the observed spectra up to 
transverse momenta of 2GeV/c, in a wider range than 
most of the previous calculations. This indicates that 
a long hydrodynamical evolution generates the correct 
amount of collective flow. For central collisions we obtain 
an excellent description of pion and kaon spectra for non- 
central rapidities. This demonstrates that the fireball is 
thermalized and that particle production is statistical for 
all rapidities in the range ±3.5 units. This confirms in 
a dynamical calculation the applicability of the statisti- 
cal fits from Ref. [77|. The elliptic flow shows a strong 
pseudorapidity dependence, that can be reproduced as- 
suming a collective thermal evolution but with reduced 
initial energy density at non-zero space-time rapidities. 

Another important result is the satisfactory descrip- 
tion of pion interferometry radii. Assuming a rapid ex- 
pansion of the system, the right amount of transverse flow 
and a reasonable life-time of the fireball are obtained. 
This gives HBT radii R ou t, Rside and Ri on g similar as 
in the experiment. The ratios Rout/Rside come out to 
within less than 10% of the measured values. We also 
show that very similar HBT radii can be obtained from 
a 2 + ID hydrodynamic calculation with correctly chosen 
initial conditions. It must be stressed however, that the 
initial energy density that reproduces the experimental 
spectra and HBT radii is 86GeV/fm 3 at the maximum 
in 2 + ID, whereas it is 107GeV/fm 3 in 3 + ID. As men- 
tioned, our 3 + ID calculation describes particle emission 
at non-zero rapidities as well. 
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